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Abstract: Truncations of the 2PI effective action are seen as a promising way of studying 
non-equilibrium dynamics in quantum field theories. We probe their applicability in the 
non-perturbative setting of topological defect formation in a symmetry-breaking phase 
transition, by comparing full classical lattice field simulations and the 2PI formulation 
for classical fields in an 0(A) symmetric scalar field theory. At next-to- leading order in 
1/A, the 2PI formalism fails to reproduce any signals of defects in the two-point function. 
This suggests that one should be careful when applying the 2PI formalism for symmetry 
breaking phase transitions. 
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1. Introduction 

An exact numerical solution of non-equilibrium dynamics in quantum field theories is only 
possible in the very simplest cases, in contrast to vacuum or thermal equilibrium. One 



generally has to resort to drastic approximations. In the classical approximation f2q , 27], 
the Hamiltonian equations of motion are solved for an ensemble of initial conditions, and 
observables are averaged over the ensemble. The approximation includes non-linear and 
non-perturbative effects, but its applicability is restricted to cases when field occupation 
numbers are large, such as high temperatures. 

An alternative approach is to study the evolution of the correlators themselves using 
Schwinger-Dyson equations in real time. These require truncation of an infinite series 
of diagrams, ordered according to some expansion of choice. One organisation of these 
diagrams follows elegantly from variations of the 2PI (or nPI) effective action, from which 
results Schwinger-Dyson equations for the 1- and 2- (and up to n-) point correlators Q. 
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Unlike naive perturbation theory, truncations of the Schwinger-Dyson equation based 
on the 2PI effective action conserve energy in time, making the approach suitable for study- 
ing the dynamics over extended time intervals (see for instance [p|). It has, for instance, 
been used to investigate equilibration and thermalisation of both scalar and fermion fields 
|3|, [|, [5], ||, [jj, reheating after inflation |8|, [| and to calculate transport coefficients [10, 11, 12 



and critical exponents [fbl. It is renormalisable [14, 15, 16, 17, 18| for all truncations, and 



is in particular applicable for very late times, where the classical approximation will break 
down as systems equilibrate classically rather than quantum Q. 

In this paper, we will concentrate on the formation of topological defects (or solitons 
fi~9f| ) in a symmetry breaking phase transition [2C, 21]. This is a highly non-perturbative 



non-equilibrium process. It also has physical relevance, because topological defects are 



produced in phase transitions of various condensed matter systems (see, for instance [21] 
and references therein) and they may have also been formed in the early universe, at the 
end of inflation p2| , |23[| . While it would be interesting to study this process in quantum 
field theory, we restrict ourselves to a classical field theory. By using the classical limit of 



the 2PI equations [ 24 ] , we can compare the results of the 2PI approach to the full numerical 
solution of the classical equations of motion, which is in principle exact. 

In a classical simulation, topological defects can be easily detected in the resulting 
field configurations, but in the 2PI approach, the basic quantity available to us is the fully 
averaged two-point function. Therefore, we first show how the presence of defects, kinks 
for N = 1 and textures for N = 2, manifests itself in the two-point function. 

In section ^ we introduce the model we are considering, an O(N) model of scalar fields 
in 1+1 dimensions. In section |3] we derive tell-tale signals in the equal-time two-point 
correlator of the presence of defects. Then in section ^| we perform detailed simulations 
of the full classical theory and establish how these signals manifest themselves. The 2PI 
formalism is then intoduced in section [| and we compare numerical simulations of the LO 
and NLO approximations to the full classical result. We conclude in section |6[ 

2. Model and defects 

We consider a classical theory of ./V real scalar fields 4> a with a G {0, ..,N — 1} in 1 + 1 
dimensions. The continuum action is 



S = Jdxdt (-8^8^ - V(<f> a )\ , (2.1) 



and has an O(N) symmetry. We shall investigate the dynamics of the system in a simple 
setup in which the potential varies with time. Initially, it corresponds to a free field with 
mass fi, 

Virata) = \^4>a4>ai (2.2) 

and at time t = 0, it changes instantaneously to 

V(<f> a ) = -\fMa + ^{Ma)\ (2.3) 
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triggering a symmetry breaking transition. We add a small constant damping term Tdt4> 
to the equation of motion, so that it becomes 

d?Mx, t) + T d t 4> a {x, t) - d 2 x <j) a {x, t) - ix 2 4> a {x, t) + ^- (Mb) (f> a {x, t) = 0. (2.4) 

biv 

This ensures that the system reaches eventually a zero-temperature state with the 0(N) 
symmetry broken spontaneously and § a §a = (6-/V/A)/i 2 = v 2 . Without damping the system 
would equilibrate in a state with a non-zero temperature, and the symmetry would be 



restored. The equation of motion (2.4) is discretised in the spatial direction on a lattice of 
spacing a and in the temporal direction using the leapfrog algorithm with time step 5t for 
the time derivative. 



For N = 1, the model has two degenerate vacua at cj>o = ±v = ±y6/A//, and there are 
topological defects, kinks, which interpolate between them. The classical kink solution is 

x 

0kinkO) = utanh-, (2.5) 

where d = \[2j \i is the kink thickness. 

For N = 2, the vacuum manifold is a circle, and the model does not have localised 
defect solutions. However, there are textures which correspond to field configurations that 
wind around the vacuum manifold. In infinite volume, they are unstable against growing 
to infinite size and becoming indistinguishable from the vacuum. However, if the system 



has a finite size L, the classical theory has stable texture solutions [25] 



*(:c) = <hix) + ifaix) = *fe(a:) = ve^ iN ^ L , (2.6) 

where the winding number N w is an integer. Because textures are stabilised by a finite 
energy barrier rather than a fundamental conservation law, they are in fact only metastable 
in the presence of thermal or quantum fluctuations. For N > 2 there are no topological 
defects. 

We choose the initial conditions at time t = to mimic the quantum vacuum state 
corresponding to the potential (|2.2|). Because this is a free theory, the equal-time quantum 
two-point functions of the field 4> a and its canonical momentum 7r a = dt4> a are simply 

(MVMq)) = (27r)5(k + q)5 ab ^-, 
Mk)ir b (q)) = (27r)5(k + q)5 ab ^, 

(Mk)n(Q)) = 0, (2.7) 



where LOk = \/k 2 + p?. Our initial conditions are given by a Gaussian ensemble of field 
configurations which has these same two-point functions. This choice of initialisation has 



been used extensively in the study of inflationary reheating in cosmology [^, 27], also 
using the 2PI formalism § In our case, it is not important how well these initial condi- 
tions reproduce the actual quantum dynamics, since we are only interested in the classical 
dynamics. For the 2PI formalism, the initial ensemble has to be Gaussian and this is a 
simple and convenient choice. 
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The classical equation of motion ( |2.4[) allows us to rescale the coupling A to unity, sug- 
gesting that only the dimensionless ratio \/fJ? plays a role. However, the initial conditions 
( |2.7| ) remove this freedom. We therefore keep both the coupling and the mass parameter. 

We solve the time evolution of the system using the initial conditions ( |2.7| ) using two 
different approaches: the full classical equations, and the 2PI formalism. In the former 
case, the classical equation of motion ( |2.4| ) is solved numerically for a large number of 
initial configurations that are drawn from the distribution specified by Eq. ( |2.7| ). Apart 
from the discretisation and finite-size errors, which should be similar in both cases, the only 
error in the full classical approach is statistical, due to a finite number of initial conditions. 
In contrast, there is no statistical error in the 2PI approach, but the truncation of the 
Schwinger-Dyson equation introduces a systematic error. 



3. Looking for defects in the propagator 
3.1 Kinks 

In classical simulations, it is easy to study topological defects by considering individual 
realisations of the field. Each realisation is inhomogeneous, and a kink, for instance, cor- 
responds to a point where <f>o vanishes. In contrast, the 2PI formalism does not give 
information about individual realisations, and if the initial ensemble is symmetric, the 
mean field remains zero {(pa) — 0. All information about the system is encoded in the 
full two-point correlator G a b(x,y,t,t') = (4> a (x,t)(j)i ) (y,t')) , which is invariant under O(N) 
transformation and spatial translations, G a b(x, y, t, t') = 5 a i,G(\x — y\, i, t'), and corresponds 
to an average over the whole ensemble. Therefore, we need to know what effect topological 
defects have on the correlator. Once we know that, we can calculate this quantity in full 
classical simulations, averaging over an ensemble of realisations, and in the 2PI formalism 
using different truncations. 



A popular approach pq , 30] is to assume that the field cp a is Gaussian. In that 



case, the density of zeros of the field is given by 



1 G"(0,M) r „^ 

Since the zeros can be identified with kinks, this would then give the number density of 
kinks, n = uq. However, since <p a becomes non-Gaussian soon after the transition, this 



approach does not work [29|. Instead, we need a method that works for strongly non- 
Gaussian fields. 



Another approach advocated in [31] is to look for a characteristic length scale L max m 
the propagator and identify it with the typical distance between kinks, so that n ~ 1/Zy max . 
It is presumably very generally true that a given number density n of kinks does indeed 
introduce some feature at the corresponding length scale. However, the existence of a 
characteristic length scale does not imply that there are kinks in the system, and this 
approach is therefore not suitable for our purposes. 

Let us, instead, calculate directly what the form of the two-point function is in the 
presence of defects. We consider a one-dimensional lattice with spacing a and assume that 
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there are randomly distributed infinitesimally thin kinks with number density n. Let us 
further assume that there are no other fluctuations, so that the field only takes values ±v. 
The changes of sign correspond to the locations of the kinks. We can choose the field to 
be positive at point x. 

Clearly G(0, t, t) = v 2 . The probability of there being a kink between points x and x+a 
is na, and therefore the neighbouring point has the value <fio(x + a) = —v with probability 
na, and 4>o(x + a) = +v with probability 1 — na. As a result G(a,t,t) = (1 — 2na)v 2 . 
At distance 2a we have <j)(x + 2a) = +v if there are either two or zero kinks between the 
points, whereby G(2a, t, t) = ((1 — na) 2 + na 2 — 2na(l — na)) v 2 . In general, 

G(Ma,t,t) = v 2 Yl (~l) fc (na) fc (l - na) M ~ k M [ = v 2 {\ - 2na) M 

k=0,..,M >' 
^0 v 2 e -2n\x-y\^ ^ 

where we have kept \x — y\ = Ma constant when taking the limit a — > 0. The Fourier 
transform of the correlator is 

An 

We therefore expect that the correlator has this form in the presence of kinks. A fit of this 
form allows us to determine the kink density n. Note, however, that since we assumed that 
the kinks are infinitesimally thin, this expression is only valid at distances longer than the 
kinks thickness, i.e., k u. 

Note also that the correlator has a power-law form G(k, t, t) oc k~ 2 at intermediate 
length scales n <ti k <C u. This is a special case of what is known as a "Porod tail" in the 
correlator. For general N and spatial dimension D, G(k,t,t) oc /^( iV + D ) ( see f or instance 

Unfortunately, Eq. ( [3.3| ) by itself is not an ideal indicator of the presence of defects, 
because the correlator of a weakly-coupled scalar field in thermal equilibrium has the same 
form, 

T 

G t herm{k,t,t) = -j- (3.4) 

Therefore, we need to take into account the finite thickness of the kinks. This gives an 
extra multiplicative factor to the correlator, 

An 

G(k,t,t) = ^ T¥ -\ ( f )kink (k)\ 2 , (3.5) 

where ^kink(fe) is the Fourier transform of the kink profile. In the absence of any fluctua- 
tions, this expression should be valid for k <ti 1/a. 
For the kink solution in Eq. ( |2.5| ), we obtain 

2iv hirkd 

0kink(AO = -r- . I i - , (3.6) 
k smli 2 vr/ca 
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and therefore we should find 



4ra / ^irkd i 2 



G(k,t,t) = — 5 pr v 2 . (3.7) 

Since d = V2//J, and t> = -y/6/A/i are known, this expression has only one free parameter n. 
In the presence of thermal or quantum fluctuations, these parameters and the kink profile 
would change. In fact, one could use Eq. ( j3lf ) to measure the kink shape </>kink(fc)- 



Using Eq. (3/?) we can also calculate, which value the Gaussian formula Eq. ( |3.1[) would 



give for the kink density. A straightforward calculation shows that in the limit nd -C 1, 



_ i fdkk*G(k) ^ 2 nr 

n °~n\j JdkG(k) ~^V3d' ( } 

Since this result is not even proportional to n but is sensitive to the kink shape, it is clear 
that Eq. (|3.1|) cannot be used to measure kink density. 



3.2 Textures 

As is obvious from Eq. (|2.6| ), textures are pure plane wave configurations, and therefore 
their signature in the two-point correlator is very simple. In a finite system with length L, 
a texture with winding N w gives a contribution Lv 2 /2 to the correlator at wave number 
k = 2irN w /L. Thus, if textures with winding N w appear with probability p(N w ) in the 
ensemble, the correlator will simply be 

G( W )=p{^)^f. (3.9) 
We can even derive the shape of the probability distribution p(N w ), if we assume that 



the textures are formed by the Kibble mechanism [20|. Immediately after the transition, 
the field 4> a is more or less constant at distances less than some "freeze-out" length scale £;, 
and completely uncorrelated at longer distances. In length L, there are therefore N% = L/£ 
uncorrelated regions, which we label by £ = 1, ... , N%. In each of them, the field has some 
constant value, 

$i = ve ie ', (3.10) 

where the phase angle 6% is random. The field interpolates smoothly between these values 
when we go from one region to the next, and it is natural to assume that it follows the 
shortest path on the vacuum manifold. The phase angle changes by an amount A#j = 
— 8i] w , where the subscript it indicates that we choose it to be in the range — n < 
A6i < it. When we follow the field throughout the whole system, through the boundary, 
back to the original point, the phase angle may wind around the vacuum manifold some 
number of times N w given by 



^A0j = 2nN w . (3.11) 
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When we let the system evolve after the transition, the winding number N w remains 
unchanged and the system equilibrates in the local minimum, which is the texture solution 
with winding N w . 

To estimate p(N w ), we note that the changes A9{ of the phase angle are random with 
a uniform probability distribution in (— w, it]. According to the central limit theorem, the 
probability distribution in the limit N% — ► oo is Gaussian, 



p < w »> = \Mr p (~^f) • <3 - 12) 

The correlator should therefore be 

G(k,t,t) = ^i/^exp ("|^ 2 ) • (3-13) 

The only free parameter in this expression is £. 

Again, we have to be careful when comparing this result with measurements, because 
one might expect the correlator to be Gaussian at very low k in any case. This is because 
massless Goldstone modes with k < T/2 are overdamped and will decay as exp(— (k 2 /T)t), 
giving 

G(k, t, t) ~ exp i-^k 2 ) . (3.14) 



r 

A Gaussian shape by itself would therefore not be evidence for the presence of textures, 
and therefore we have to check that the exponent approaches a constant at t — ► oo rather 
than growing linearly as Eq. (3.14) would predict. 



4. Full classical simulations 

We discretised the equations on a lattice with spacing a = 1, which means that all dimen- 
sionful quantities are expressed in lattice units. The time step was St = 0.1. In all runs, 
we used the coupling A = 0.6. 

4.1 Kinks 

To test the prediction ( |3.7D we carried out 2000 runs with ii 2 = 0.01 on a lattice with L = 
16384. The damping rate was T = 0.001. The kink thickness in this case is d ~ 14.14 ^> 1, 
and the kinks should therefore be well approximated by the continuum solution fl2.5[ ). 
In Fig. m we show the two-point function measured at time t = 10000, together with 



a fit of the form (3.7). The only free parameter in the fit is the kink number density 
n = 0.00116(1). Only points at k < 0.15 were included in the fit. At higher k, the correlator 
is still dominated by the perturbative "quantum" fluctuations, which decay exponentially 
with time because of the damping term. 

The fit is very good at k < 0.15 (i.e., k/fi < 1.5), which shows that the correlator is 



dominated by kinks of the form (2J3). Further evidence for this is shown in the inset of 
Fig. m, where we compare the number density of kinks determined by fitting the measured 
two-point function with Eq. ([O]) to the result obtained by counting the zeros of the field 
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Figure 1: The two-point function G(k, t, t) for N = 1 at time t = 10000. The dashed curve is a fit 

116(1) x 10~ 5 . The inset shows the number density of kinks measured 
directly from field configurations (solid line), fitting 



of the form (3.7) with n 
in three different ways: Counting zeros of 



Eq. (3.7) to the measured correlation function (dashed line) and using the Gaussian formula in 



Eq. fly. 



(dotted line). While Eq. ( p.l| ) works well at early times when the field configuration is 



Gaussian, it fails at later times when the kinks have actually formed. In contrast, Eq. (3.7) agrees 
very well with the direct measurement at late times. 



in the lattice field configurations. The two results agree perfectly. For comparison, we also 
show the incorrect result obtained with Eq. fl3.ip , 

4.2 Textures 

To test Eq. (pp , we chose L = 4096, /j 2 = 0.49 and T = 0.01. Fig. ^ shows the correlator 
at various times, together with a fit of the form ( |3.13| ). Again, the fit is very good. The 
inset shows the fit parameter £ at various times, together with an exponential fit 



C(t) = £(oc) - A£e 



-at 



(4.1) 



The fit parameters are £(oo) = 21.4(2), A£ = 18.4(1) and a = 2.68(6) x 10~ 4 . This shows 
that the length scale £ approach asymptotically a constant at late times, ruling out the 
diffusive behaviour of Eq. (|3.14 ) and providing a clear signal for the presence of textures. 



4.3 Higher N 

For N > 2, there are no topological defects, and the system should end up in the vacuum 
state, in which the two-point function is simply a delta function. In finite volume, 



lim G(k,t,t) 



v 2 L/N for k = 0, 
fot A; j=- 0. 



(4.2) 
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Figure 2: The two-point function G(k,t,t) for N = 2 at times t = 1000 (white), t = 3000 (grey) 



and t = 10000 (black). The curves are Gaussian fits of the form fl3.13| ) with £ = 7.23 ± 0.02, 
£ = 13.03 ± 0.04 and £ = 19.94 ± 0.22. In the inset, the data points show the time evolution of the 
fit parameter £, and the curve is an exponential fit. 



In Fig. |, we show the correlator for N = 4. It appears to be well fitted by an exponential, 

G(k,t,t) oc e~ lk . (4.3) 

As the inset shows, the exponent grows as I oc i 2 / 3 , becoming singular at t — > oo as 
expected. 



5. 2PI-1/N approximation at LO and NLO 

The 2PI formalism is based on quantum field theory, but it can be used to study an 
ensemble of classical field configurations. The corresponding equations are obtained by 
taking the classical limit h — > of the quantum equations. In our case, this is not entirely 



straightforward because our classical equation of motion (2.4) contains a damping term, 
and it is therefore not immediately obvious what the corresponding quantum theory should 
be. 



To avoid this problem, we note that Eq. (2.4) is also the equation of motion of an 
undamped 2+1-dimensional scalar field theory in an expanding anisotropic space with 
metric 



ds 2 



dt 2 - dx 2 - a(t) 2 dy 2 



The equation of motion for the scalar fields <p a is 



dfya + -d t <Pa - 8^(f> a (x,t) - \dl(f) a (x,t) + V'{4> a ) 



(5.1) 



(5.2) 
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Figure 3: The two-point function G(k,t,t) for N = 4 at times t = 1000 (white), t = 3000 (grey) 
and t — 10000 (black). The curves are exponential fits to the first six data points. The plateau (or 
k^ 1 behaviour to be more precise) at higher k is a remnant of the initial "quantum" fluctuations 
(2.7). In the inset, the data points show the time evolution of the exponent, and the blue curve is 
a power-law fit with exponent 0.670(6). 



If the scale factor grows exponentially, 

a(t) = a(0)e Tt , (5.3) 

and all the fields 4> a are independent of y, this reduces to Eq. (2.4). However, since the 
dynamics is undamped, the quantum generalisation of the system is obvious. 
To derive the 2PI equations, we write the action as 1 

S = J a(t)dydxdt (^d t (f> a d t (t> a - ^d x <j) a d x <f) a - 2 ^(t) dy ^ ady( ^ a ~ V ^ a ^j ' ( 5 ' 4 ) 

Imposing that the field is explicitly independent of the y-coordinate (f>(x, y, t) = cp(x, t), we 
have 

S = L, J a^it (I 8A8A - - r<u) , (5.5) 

where we have performed the integration over y, J dy = L y . We end up with an effective 
theory of N real scalar fields in 1 + 1 dimensions, but with a non-standard action, including 
a time-dependent factor a(t). In the classical limit, the factor L y drops out, so we can ignore 
it. 



1 The action is in principle defined on the Keldysh contour S = Sc- We will suppress this complication 
for the moment, since it has no implications when deriving the classical equations of motion. It does enter 
and is crucial in the derivation of the 2PI evolution equations (see also Appendix |B|). 
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We note that we can write the action as 

S = J dxdt ^-a{t)<P a [d 2 t +Td t -dl-^]^ a -^{^ a ?y (5.6) 

In the 2PI formalism, 2 evolution equations for the correlator G(x,y,t,t') are derived 
by varying the 2-point irreducible (2PI) effective action corresponding to Eq. (|5.5| ). The 
result is the Schwinger-Dyson equation for the statistical (F) and spectral (p) parts of the 
propagator G(x,y,t,t') = F(x,y,t,t') — |sign(t, t')p(x,y, t, We assume that the mean 
field is zero and explicitly impose 0(iV)-symmetry G a j> = 5 a t,G, as well as homogeneity and 
isotropy, G(x, y, t, t') = G(r, t, t'), r = \x — y\, so that 

(d 2 - d 2 x + Td t + M 2 (t)) F(r,t,t') = - [ dz dt" S p (r - z,t,t")a(t")F(z,t" ,t') 

Jo 

+ / dzdfZptr - z,t,t")a(t")p(z,t",t'), (5.7) 
Jo 

(dt -dl + Td t + M 2 (t)) p(r,t,t') = - [ dz dt" S p (r - z,t,t")a(t")p(z,t" ,t'), (5.8) 

Jt' 



3N I' 1 

—D F (r,t,t') = -TL F (r,t,t') + / dz dt" U p (r - z,t,t")a(t")D F (z,t" ,t') 
* Jo 

- / dzdt"U F (r - z,t,t")a(t")D p (z,t",t'), (5.9) 
Jo 

—D p (r,t,t') =-U p (r,t,t') + / dzdt"U p (r- z,t,t")a(t")D p (z,t",t'), (5.10) 
^ Jo 



with the effective mass given by the local part of the self-energy 



M 2 (t) = -p 2 + A^±^F(0,M). (5.H) 

DiV 



The objects Dp and D p are components of an auxiliary propagator [33]. The non-local 
parts of the self-energies are given at NLO of the 1/N expansion as 

EHr, = -3^ [Hr, t, t')D F (r, t, t')] , (5.12) 

S p (r, t, t') = - A [p(r, t, t')D F (r, t, t') + F(r, t, t')D p (r, t, t')) , (5.13) 

n F (r, t, t') = -j [F(r, t, t')F(r, t, t')] , (5.14) 
U p (r, t, t') = —NF(r, t, t')p(r, t, t'). (5.15) 



We have taken the classical limit, as prescribed in J24} |. 

The functions F(p), Y, F (T, p ), H F (H p ) and D F (D p ) are (anti-) symmetric in t,t' and in 
particular p(t, t')\ t=t > = 0. When T = we should set 

d t p(tj)\t=t> = 1, (5.16) 

2 For details of the implementation of the 1/N expansion applied to the O(N) model see H |3^] and 
Appendix H. 
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for all times. With the damping term, we find instead 



d t p{t,t% =t , oce- r *. (5.17) 

The scale factor a(t) only appears in the self-energy, for every time slice t" taking on the 
appropriate value a(t"). 

At leading order (LO), N — * oo, we should drop all the non-local self-energies (right- 
hand sides of Eq. (|5,7| )) and take the limit 

N + 2 1 

-W - 5- (5 ' 18) 

in Eq. ( |5.11 ). Keeping all of ( 5.11 ) constitutes the Hartree aproximation, which is leading 



order in a coupling expansion, but a mixture of LO and some NLO in the 1/N expansion. 
Clearly, since at LO the equations have no reference to N, it is impossible for the approx- 
imation to know about defects. The Hartree approximation amounts to a rescaling of the 
coupling A — > X(N + 2)/N, and so is equivalent to LO 3 . 

5.1 Spinodal transition in the free and LO approximations 

At very early times, we can neglect the coupling altogether, and solve for the correlator as 
the spinodal transition proceeds, 



-T(t+t')/2 



4< 



ib 2 \ ( U! 2 \ 

1 + -§ cos[cj fc (i - t')} + 1 - -|) cos[uj k (t + t')} 

-2% ( sinL>fc(£ - 1')] + it sm[uj k (t + 1')] J . (5.19) 
\u k )_ 

where (see Appendix A) , uoj: = \/k 2 rb /i 2 , and 

L0 2 k = k 2 -{ f i 2 + r 2 /A), ^ = fc 2 +( / x 2 + r 2 /4), f = -^. (5.20) 

Fig. ^ (left) shows the evolution of the equal time, equal space correlator at early times 
in the free, LO and NLO approximations, as well as for the full classical simulation. The 
initial growth is faster than exponential. The back-reaction of the interaction term ends 
the growth at (N)LO, and the correlator starts oscillating around it's vev. The apparent 
damping at LO is the effect of different frequency modes coming out of phase ( "dephasing" ) . 
No memory is lost through dissipation, and (partial) recurrence is seen at later times, when 
the (most dominant) modes come back into phase. At NLO the damping is real, and the 
system will eventually thermalise 

Fig. |] (right) shows the equal time correlator in momentum space as the tachyonic 
transition takes place. Around time t = 4, the back reaction kicks in and growth ends. 
There is a clear separation between unstable k < fi and oscillating modes k > fi. 



3 Note, that in the presence of a mean field {(j>) — v, at LO the fluctuations around this mean field 
satisfy Goldstone's theorem (and so has zero modes for N > 1), whereas the Hartree approximation does 
not. When (</>) = 0, Hartree acts as LO with the rescaled coupling, and hence finds a different value for 
v 2 — > jf^v 2 (F(r — 0) — v 2 N/(N + 2)) around which it again has zero modes. 
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Figure 4: Left: Equal time correlator F(t,t,r = 0) in time for the free, LO, NLO approximations 
and the full classical simulation. Right: The correlator \n(G(k,t, t)) vs k in the free (solid lines) 
and LO (dashed lines) approximations. n x — 256, A = 0.6, T = 0.02, [x 2 = 0.49. 



5.2 Numerics at NLO 



Whereas the LO correlator only knows about the massless (Goldstone) modes, at NLO the 
correlator holds information of both the massive (Higgs) mode and the massless modes. 
The correlator turns out to be the average of a massive correlator and N — 1 massless 
correlators (at late times) ||. In this sense, NLO is a great improvement on LO. 

We carried out simulations for N = 1,2,4 and 16 using both the full classical equations 
and the NLO 2PI approach. The parameters were identical in both sets of simulations: 
A = 0.6, n 2 = 0.49, T = 0.02, 5t = 0.1, and the lattice size was L = 256 (in lattice units). 
In the classical simulations, we averaged over 2000 different initial conditions. This was 
enough to reduce the statistical error to such a low level that it is impossible to see the 
error bars in any of our plots. 

Fig. [5] shows the equal-time momentum-space two-point function G(k) at an early 
time t = 10. Apart from a small difference at low k in the N = 1 data, the 2PI approach 
reproduces the classical results very well. 

In Fig. ||, we show the corresponding plots at a later time t = 190, when the system is 
approaching the equilibrium state. While N = 4 and N = 16 still agree very well, there is 
a clear discrepancy between the NLO 2PI and the full result for N = 1 and N = 2. 



For iV = 1 we did a one-parameter fit to the full data using Eq. (3.7) with the kink 
density n as the only free parameter. As the plot shows, this fits the data very well with 
n = 0.016. Therefore, we have to conclude that the contribution from the kinks is missing 
from the NLO 2PI result. 



Similarly, we fitted Eq. ( 3.13 ) to the N = 2 data using the length scale £ as the only 
free parameter, and obtained a very good fit with £ = 13.0. The hump in the correlator at 
low k is therefore due to textures. Again, there is no sign of this contribution in the NLO 
2PI data. 
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Figure 5: The equal-time momentum-space correlator G(k) at an early time t = 10 in the full 
classical simulation (dots) and the 2PI-NLO approximation (line) for different N. The parameters 
were A = 0.6, = 0.49, T — 0.02. The agreement is very good for N > 1 and reasonable even for 
N = 1. 



6. Conclusions 

We have shown that the 2PI formalism at next-to-leading order in a 1/N expansions fails to 
describe topological defects, which exist in the 0(N) scalar field theory when N < 2. This is 
disappointing, because this method has been seen as a promising way of addressing generic 
non-equilibrium questions in quantum field theory. Instead, one has to be very careful 
when using the method and make sure that all relevant effects are included, especially 
since there is no hint in the NLO 2PI results themselves that they break down. 

One possible caveat is that we are applying a 1/N expansion at low N. Still, the 
fact that the discrepancy is qualitative at both N = 1 and N = 2, specifically in the 
momentum range sensitive to defects, and small already at N = 4, leads us to believe that 
the shortcomings are not a result of the choice of expansion. Going to higher orders in the 



1/N approximation 34] is not likely to improve the situation, because each type of defects 
is specific to one particular value of N, and their contribution is therefore non-analytic 
in 1/N. Unfortunately, the obvious alternative, the loop expansion, is in fact numerically 
unstable in the presence of large occupation numbers as in the present case of a spinodal 
transition. 

What this means for the 2PI formalism is that as for any perturbative framework, 
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Figure 6: The equal-time momentum-space correlator G(k) at a late time t — 190 in the full 
classical simulation (dots) and the 2PI-NLO approximation (line) for different N. The parameters 
were A = 0.6, [x 2 = 0.49, T = 0.02. The agreement is very good for N = 4 and N = 16, but 
for N = 1 and N — 2 there is a clear discrepancy, which matches precisely the predicted defect 
contribution (dashed lines). 



certain observables are beyond the reach of the approximation. One should be careful when 
studying (symmetry breaking) phase transitions, and be aware that quantities that are 
sensitive to the presence of defects will come out wrong. Indeed, since defects (monopoles, 
vortices, strings, textures, instantons) contribute to the propagator even in equilibrium, 
certain phenomena may not be reproduced by the 2PI formalism. Since the expansions on 
which current methods are based are unlikely to describe topological defects correctly, one 
may need a radically different approach, perhaps along the lines of the Hartree ensemble 
method developed in Ref. ]35| ], 
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A. Free field 

It is useful to take a step back and take a look at the simplest approximation to the 
dynamics, the free theory. This corresponds to very early times, as the potential is quenched 
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from positive to negative curvature. The field experiences a spinodal instability, and we 
can solve exactly for the field evolution and the correlator in the potential 

f 



V(4> a ) = V - -^4>a4>a- 



(A.l) 



The classical and quantum equations coincide, and we can think in terms of the free Klein- 
Gordon equation, which reads 



(d 2 t +rdt + ^> 2 ) k = o, 



with uj^ = \/k 2 ± /i 2 . The solution is of the form 

Mt) = e~ Tt/2 (oke^ 4 + /? k e-^) , 

with the momentum 

7Tk(t) = d t cf> k (t) = ioJke-'^ 12 

We use 



1 + if) a k e iuJkt - (l - iT\ (3 k e 



l = k 2 -{^ + r 2 /4) , 



~.2 



k 2 + (n 2 + r 2 /4) , f 



2uj k 



(A.2) 

(A.3) 
(A.4) 
(A.5) 



Before the quench at t < 0, the system is in the vacuum state in the potential V 
b 



Vq + hn 2 4> a (p a , so that 



0k (o- 



1 icu^~ 

^(a k + al k ), 7r k (0~) = r L_ (a k - al k ), 



2w ( 



(A.6) 



2^7 



with (a k aj) = <5 3 (k — 1), and all other correlators of a's and a^'s are zero. 

Matching Eqs. ( |A.3| ) and (A.4) to (A.6) at i = 0, it is then straightforward to calculate 
the correlators of interest, 



-T(t+t')/2 
4w+ 



1 + ^| ) cos[w fc (t - t')] + ( 1 - ^| ) cos[w fc (t + t')] 



^ 2 



-2i ( ^- smL> fc (i - i')] + »f sin[o; fc (i + t')} 



^(i)4(t')-7T k (t')4(*)} 



-ze 



-r(i+i')/2 
-r(t+t')/2 



cos[a;fc(t — t')] + f sin[a;fc(t — t') 



(A.7) 
(A.8) 



2<4 



Ll)J 



r 1 + i (cos[co k (t - t')] - cos[u k (t + t')}) 



Co 2 



+ 1 + -| sm[to k (t - t')] - 1 - -§ + 2f 2 sinK(t + i' 



(A.9) 



M*)4(0> 



2 _-r(t+t')/2 



4<4 



2iu; 



+ 



r 2 + 1 1 + -| cos[io k (t - t')] + sin[w fc (t - t')] 



u k 



,~,2 



LO'7 



+2 r T 2 --| sin[u; fc (i + 0]+ r 2 -i 1 



4f 2 J cos[w fc (t + t')] 



(A.10) 
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It is easy to convince oneself that these correlators match the results of in the limit 
T = 0,t = t'. 

Modes with k > \f n 2 + T 2 /4 simply oscillate in time. For modes with k < yV 2 + T 2 /4, 
the oscillation turns into exponential growth (for large t ±t'), 

(0k4)(i) oc e~ v ^' 2 ( e KK*+*') ± e KK*-*')) . (A.ll) 

For equal time t = t', we have 

(0 k 4)(t) oc exp (2t Wpfl + ra/4 — k 2 — T/2]) . (A.12) 



So for all r, modes with k < fi (rather than k < 1//X 2 + T 2 /4) are unstable. Also, 



V 1+ (^) "i^^ 1, r ^°' ^4)Woce 2 KTl*, (A.13) 

— > 0, T — > oo, (0k<^k)(i) « constant. (A. 14) 

It is of course straightforward to generalise the initial state to a state with arbitrary quasi- 
particle numbers n^, 

(a k a k ) = n k + 1, (a k a k ) = n k , (A.15) 
for instance starting from a thermal initial state 

n k = (e u ' k l T — 1^ , quantum, n k = T/uj^ classical. (A. 16) 

B. 2PI equations of motion 

The 2PI equations result from variation of the 2PI effective action, which reads |l|] (sup- 
pressing space labels) 

T[G] = l -Tr In G- 1 + '-G^iG - Go) +T 2 , (B.l) 

where 

and Sq is the free action. iT2 is the sum of vacuum diagrams in terms of the Feynman 
rules, starting at two loops. Extremisation gives the physical propagator 

JIt -™^ 1 - 2 ^- (B - 3) 

Multiplying from the right by G, we have 

J dt" (Go 1 (t,t")-E(t,t"))G(t",t') = 8(t-t'), (B.4) 
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where all integrals are to be performed along the Keldysh contour C. In the case of our 



action Eq. (5,5) 



G 1 (t , t') = ia{t) (d t 2 - 8 2 x + T d t + /x 2 ) . (B.5) 
The coupling is effectively time-dependent \{t) = Xa(t), so that the self-energy is given by 

5 2iT 2 [G,\(t)} 

S(M) = 6G(t,t>) ' (R6) 

For the purpose of illustration, we will look at the coupling expansion. At order A, the 
effective action contains a single diagram, the "figure-8", 

iT 8 2 = -i\(t) J^dtG 2 ab (t,t). (B.7) 

This leads to the Hartree approximation, for which we find 

a(t) (a 2 - d 2 x + Td t - fi 2 ) G(t, = -X( t )^^F(t, t)G(t, t'). (B.8) 

In this case, a(t) cancels out, and the equation of motion reduces to the usual Hartree 
approximation, with an added damping term. This makes perfect sense, since at Hartree 
order, classical and quantum evolution coincides, and a(t) appears only in the classical 
equation of motion Eq. ( |2.4| ) through the damping term. 

Beyond the Hartree order, the modification by including the expansion is to make the 
replacement A — > \a(t) for all vertices of the self-energy. So, for the sunset diagram (order 
A 2 (t)) we should add 

■psunset = (r£? [ dtdt'a{t)G\t,t')a{t'), (B.9) 
48 Jc 

£ sunsct (M') = -±- a (t)G{t,t'fa(t'). (B.10) 
6 

In the equation of motion, we add to the right hand side of Eq. (|B.8| ) 

J dt"E(t,t")G(t",t') = J dt"G{t,t"fa{t")G{t"J). (B.ll) 

Note, that a factor of a(t) will again cancel out, leaving only the a(t") which is integrated 
over. 

An alternative interpretation of this is to return to the original 2+1 dimensional sys- 
tem Eq. (|5.4|). If we consider the scale factor a(t) as part of the metric rather than the 
Lagrangian, then A is time-independent, and the integrals become 

J dt" dxa(t")dyX 2 f(t,t',t",x,y) = J dt" dxL y a(t)X 2 f(t,t' ,t" , x), (B.12) 

so that the occurrence of a(t") is a result of integrating over a homogeneous slice in y, with 
increasing size L y a(t"). 
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We can attempt a naive a(i)-power counting in perturbation theory. The free propa- 
gator Eq. ( |A.7|) is of the form 



G(t,t',a(t),a(t')) = G ffl oc -^J^, (B.13) 
^Jaifyaif) a(t) ]j a(t') 

neglecting the instability of some of the modes, for which the growth is determined by 
k 2 — fi 2 and not a(t). 

We can now write the Schwinger-Dyson equation in powers of a(t) by assuming that the 
full propagator supplies the same power of a(t) as the free G(t,t') oc a~ l {f){a(t) / a(t')) 1 / 2 , 

(0[Free] + 0[Hartree] + 0[NLO]) G(t, t') = 5{t - t'), 

(B.14) 

(OHt^+O^it)}) G(t,t') + J dt"0[a- 1 (t)^^]G(t",t') = 6(t-t'). 

(B.15) 

For each extra vertex, we get another order a -1 . This is because each extra vertex 
supplies one factor of A(i'") and four extra propagators legs (for a four- vertex) , giving 
(a _1 / 2 (t w )) 4 a(t /// ). In this way, each subsequent order of A gives an extra integration and 
a factor of a -1 . In contrast, had we simply counted orders of X(t) without reference to the 
suppression of the propagators, we would have had a factor of a(t) for each extra vertex, 
which in this case is exponentially growing in time, making the expansion unreliable at 
best. 

In practice, it turns out that the 1/N expansion is more robust towards very large 
occupation numbers, as in the present case of a spinodal instability. The inclusion of the 
expansion in the 1/N case proceeds along the same lines as for the coupling expansion 
discussed here, resulting in Eqs. fl5Y7 - |5,15| ), 
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